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Abstract. We study the fidelity decay of the A:-body embedded ensembles of random matrices 
for bosons distributed over two single-particle states. Fidelity is defined in terms of a reference 
Hamiltonian, which is a purely diagonal matrix consisting of a fixed one-body term and includes 
the diagonal of the perturbing fc-body embedded ensemble matrix, and the perturbed Hamiltonian 
which includes the residual off-diagonal elements of the A:-body interaction. This choice mimics the 
typical mean-field basis used in many calculations. We study separately the cases k = 2 and 3. We 
compute the ensemble-averaged fidelity decay as well as the fidelity of typical members with respect 
to an initial random state. Average fidelity displays a revival at the Heisenberg time, t = tn = 1, 
and a freeze in the fidelity decay, during which periodic revivals of period tn are observed. We 
obtain the relevant scaling properties with respect to the number of bosons and the strength of the 
perturbation. For certain members of the ensemble, we find that the period of the revivals during the 
freeze of fidelity occurs at fractional times of f#. These fractional periodic revivals are related to the 
dominance of specific £-body terms in the perturbation. 
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1. INTRODUCTION 

Fidelity, also named quantum Loschmidt echo, is a measure of the sensitivity of the 
dynamics of quantum systems to perturbations and has attracted a lot of attention in the 
last years in connection with quantum information processes and quantum chaology; for 
a recent review see [1]. 

Fidelity compares the time evolution of a given initial state under a reference Hamilto- 
nian with the time evolution of the same state under a slightly different one. We denote 
as Hq the reference or unperturbed Hamiltonian, and as H £ = Hq + eV the perturbed 
one which includes the residual interaction. The former is associated with the unitary 
time-evolution %j(t) = fexp[— ^ /q cItHo(t)], where f is the time-ordering operator; 
likewise, % £ (t) is the propagator associated with the perturbed Hamiltonian. In the defi- 
nitions of the propagators we have considered the most general case in which the Hamil- 
tonians may depend on t explicitly. With these definitions, and taking an arbitrary initial 
state l^o), the fidelity amplitude is defined as 

/ e (0 = <*o|^o(-0^(OI*o> s (1) 
whose square modulus is known as the fidelity 

F £ (t) = \f £ (t)\ 2 . (2) 



Fidelity decay can occur in many different ways and in successive combinations of 
these. For an extremely short time, the so-called Zeno time, all systems will decay 
quadratically in time, both for chaotic and integrable systems. After the Zeno time 
systems specific correlations often dominate the decay. 

For regular systems the decay typically changes into a quadratic and later Gaussian 
decay dominated by the dimension of the effective Hilbert space that is accessible to 
the wave function. Furthermore revivals are very common, particularly if the system has 
few degrees of freedom. 

For chaotic systems on the other hand we expect, after the system specific behavior, 
linear or exponential decay of the fidelity which is determined by the coupling strength. 
This phase is often known as Fermi Golden Rule decay. Around the Heisenberg time the 
exponential decay changes over into a quadratic function or a Gaussian. An alternative 
behavior, known as Lyapunov decay, can occur with rather strong perturbations. For 
times of the order of the Ehrenfest time, exponential decay can be observed. The 
Lyapunov exponent, and not the perturbation strength, define the decay rate. 

An alternative and interesting behavior is found if the diagonal part of the perturbation 
(or at least its time-average) vanish in the basis in which the unperturbed Hamiltonian 
is diagonal. Then we get a fidelity freeze, which has been shown for classically inte- 
grable [3] and chaotic [4] systems, as well as for random matrix models [5]. Both, the 
value at which the freeze occurs and the length of its duration, are functions of per- 
turbation strength. For fermions, it was proposed that a mean field Hamiltonian as the 
reference system and the residual interaction could well lead to a fidelity freeze [6]. 
Specifically, this problem was treated in the framework of an embedded two-body ran- 
dom ensemble for fermions [7]. While there is no freeze on average, because of long 
tails in the distribution, typical samples from the ensemble display the freeze and indeed 
the median does [6]. 

In this paper we study the fidelity decay for the fc-body embedded ensembles of 
random matrices for bosons, which are distributed over two single-particle states. We 
shall consider a reference Hamiltonian which consists purely of a diagonal 1-body 
terms as well as the diagonal fc-body matrix elements, with k = 2 or k = 3, which 
are the physical important cases. The perturbed Hamiltonian will include a traceless 
residual interaction, consisting of the off-diagonal matrix elements of the same fc-body 
interaction. Our results indicate that ensemble-averaged fidelity decays quadratically on 
time for short times (before the Heisenberg time t H ) as well as for very long times, 
precisely after the freeze of fidelity ends. Furthermore, freeze of fidelity is observed for 
the ensemble-averaged fidelity, and during the freeze, it displays periodic revivals with 
the periodicity given by the Heisenberg time. These results occur often but not always, 
when individual members of the ensemble are considered. We present numerical results 
showing that the periodicity of the revivals is an integer fraction of tn, and demonstrate 
that the specific value of this period is related with specific off-diagonal terms of the 
fc-body interaction considered in the perturbation. 



2. THE fc-BODY TWO-LEVEL BOSONIC ENSEMBLE OF 

RANDOM MATRICES 



The fc-body Embedded Ensemble of Random Matrices for bosons considers all possible 
fc-body interactions among spin-less n-boson states, where the bosons are distributed 
over / single-particle states [8]. Below we discuss the simplest case / = 2. To define 
this ensemble, we first introduce the single-particle states associated with the operators 
a| and dj, with j = 1 , 2 (7 = 2), which, respectively, create or annihilate one boson 
on the single-particle level j. These operators satisfy the usual commutation relations 
for bosons, i.e., [a;, at] = <5,j. The normalized n-boson states are denoted as |jU r ^) = 

{jV r ^)~ l {d\) r {d1 2 ) n ~ r \Qi), where jY^ = [r!(n — r)!] 1 / 2 is a normalization constant and 
|0) is the vacuum state. The Hilbert-space dimension is thus N = n+l. Then, in second- 
quantized form, the most general fc-body interaction of n bosons distributed over two 

single-particle levels, Hjf', can be written as [9] 



MP) v V W (*?) r (ffi- Wte)*"' m 



t(p) _ 



Here, /: denotes the rank of the interaction, 1 < k < n, and are the fc-body matrix 
elements, which are independent Gaussian-distributed random numbers with zero mean 
and constant (fixed) variance Vq = 1. As in the case of the canonical random matrix 
ensembles (see [10]), Dyson's parameter /3 distinguishes the cases according to time- 
reversal invariance: jS = 1 is the case where time-reversal invariance holds, and /3 = 2 
where this symmetry is broken. Hence, the £-body interaction matrix v^) is a member 
of the Gaussian orthogonal ensemble (GOE) for /3 = 1, or Gaussian unitary ensemble 

(GUE) for j8 = 2. Note that ' commutes with the number operator h = a[a\ -\-a 2 ci2, 
i.e., the interaction preserves the total number of bosons n. 

This ensemble presents some noteworthy properties: It exhibits non-ergodic level 
statistics, i.e., spectral or ensemble unfolding do not yield the same results [8]. In 
addition to this, for /3 = 1 the ensemble displays a large and robust quasi-degenerate 
portion of the spectrum for a wide interval of k, while for /3 = 2 only seldom accidental 
quasi-degeneracies are observed [11]. These results are due to the fact that each member 
of the ensemble is Liouville integrable in the semiclassical limit, for all values of k; see 
Ref. [12] for details. 



3. NUMERICAL RESULTS ON FIDELITY DECAY 

3.1. Definitions 

In order to study fidelity, as stated above, we must first define the reference (or 
unperturbed) and the perturbed Hamiltonians considered. The reference Hamiltonian 
J#o k shall be defined as the sum of a fixed one-body interaction term, and the diagonal 
part of the fc-body interaction, with either k = 2 or k = 3. Since each of these terms 




FIGURE 1. Comparisson between the time evolution of the ensemble-averaged fidelity illustrated with 
1 — (F(f )) when the diagonal part of the £-body interaction, which is included in the reference Hamiltonian 
J%o k , is fixed (upper curve for long times) or varied with the realizations (lower curve for long times) for 
different realizations of the ensemble. 



have different spectral widths [8], we normalize each one with the respective width 
of the spectrum Wk- This mean-field unperturbed Hamiltonian is further restricted by 
imposing that J%k,o is diagonal in the occupation number basis, which we shall use as 
our reference basis. Notice that the reference Hamiltonian contains the diagonal part 

of the /c-body interaction, i.e., it includes the /c-body matrix elements of the form v 
Then, the unperturbed Hamiltonian is explicitly given by 



r.r ■ 



Jtc\ 



0/, 



1 ^diag A ^-diag 

w k=l k = l w k k 



(4) 



In Eq. (4), the width of the spectrum of the /c-body embedded ensemble is given by 



N 



tr[#f))]2 



N 



(5) 



where 



n — s 



n + s+l 
k 



(6) 



In these expressions we have used explicitly the restriction to the two single-particle 
level case (Z = 2); the label B stands for bosons, the over-line indicates ensemble average, 

(s) 

and A y B '(k) is the 5-th eigenvalue of the ensemble-averaged correlation matrix of the 
bosonic /c-body embedded ensemble; see Ref. [8] for details. 

As for the residual interaction, it simply consists of the remaining off-diagonal matrix 
elements of the /c-body interaction, properly normalized by Wk. Denoting by X the 



perturbation strength, the perturbed Hamiltonian is written as 

A k = M k + ^ m ^. (7) 

Notice that with this definition of the reference Hamiltonian and the traceless residual 
interaction, where J#o k also depends upon X, the conditions to observe fidelity freeze 
are fulfilled [5]. 



3.2. Ensemble-averaged fidelity decay 

Here, we present numerical results for the ensemble-averaged fidelity decay in terms 
of the perturbation strength X and the number of particles n. For the numerical treatment, 
we have fixed the diagonal part of the fc-body perturbation to a specific realization of 
the ensemble. Therefore, the different realizations of the ensemble will only involve 
different realizations of the off-diagonal part of the perturbation. This was done for 
numerical convenience. We emphasize that this implementation has no effect in the 
results; this is shown in Fig 1, where the two curves included are difficult to distinguish 
quantitatively. In order to obtain the physically relevant scalings, it is convenient to 
rescale the physical time t' by a dimensionless time measured in units of the Heisenberg 
time, i.e., t = t'/tn- Here, tn = 2nh/d and d denotes the average level-spacing of the 
spectrum of the unperturbed Hamiltonian J#o k . 

The results for ensemble average were obtained as follows: We first fixed the unper- 
turbed Hamiltonian J#o k , and considered 1000 independent realizations of the perturba- 
tion (involving only the off-diagonal part of the fc-body interaction); for a given perturba- 
tion strength X and fixing the number of bosons n, this defined the perturbed Hamiltonian 
Mx k through Eq. (7). For each realization of the residual interaction H k , we calcu- 
lated the corresponding time evolutions from an initial random state, which permitted to 
obtain the fidelity for this particular realization, using Eq. (2). Averaging the resulting 
fidelities over different realizations constitutes the ensemble-averaged fidelity. 

In Fig. 2 we present numerical results for the 1 — (F(t)} as a function of time 
(in Heisenberg time units), as convenient representation of ensemble-averaged fidelity 
decay, for various values of the perturbation strength X (for n = 1000), and various 
values of the number of particles n (for X = 10~ 6 ). The cases illustrated corresponds 
to k = 2 and /3 = 1 ; similar results were obtained for k = 3 as well as when the time- 
reversal symmetry does not hold (/3 = 2). As shown in Fig. 2, for times smaller than the 
Heisenberg time t <tfj = 1, but longer than Ehrenfest time, fidelity displays a quadratic 
decay in time, which is the typical situation observed for integrable systems. At the 
Heisenberg time (t = tn = 1), the system exhibits a revival, i.e., {F(t)) approaches the 
unity again. This revival is natural for integrable Id systems [1], which corresponds to 
this case, but is also observed for the Gaussian ensembles of Random Matrices [13], or 
in models with underlying chaotic dynamics [14]. 

Beyond the Heisenberg time, fidelity remains essentially constant at certain value 
Fplateau'i this is the fidelity freeze. The value F p i ateau scales as A 2 and n 1 with respect to 
the perturbation strength and the number of particles, respectively. We note that during 
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FIGURE 2. Time evolution of 1 — (F(t)) for k = 2, j3 = 1. (a) Dependence upon the strength of the 
perturbation (n = 1000 particles); (b) dependence upon the number of particles for X = 1 X 10~ 6 . For 
very short and very long times, the fidelity (F(t)) scales as t 2 . (c) Illustration of the periodic revivals 
displayed by the ensemble-averaged fidelity during the freeze. Note that the periodicity is the Heisenberg 
time. 



the freeze, fidelity displays some very short periodic revivals, whose period is precisely 
the Heisenberg time tn = 1 • To the best of our knowledge, such periodic revivals during 
the freeze of fidelity have not been observed. The freeze of fidelity lasts until some 
ending time t e is reached, which displays a A -1 dependence. After t e , fidelity decays 
once again quadratically in time, until it saturates near F = 0. 
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FIGURE 2. (Continued) 



3.3. Fidelity decay of specific members of the ensemble 

Fidelity, for most of the realizations of the ensemble, resembles qualitatively the 
ensemble-average fidelity. Yet, for some specific realizations the results do differ. A 
common case is that the position of the plateau may be slightly shifted up or down with 
respect to the average one. 

A more interesting case occurs for some specific realizations of the ensemble, where 
the differences are more subtle and far reaching; see Figs. 3 and 4. Indeed, for some 
specific members of the ensemble, fidelity displays periodic revivals during the freeze 
with an integer fractional period. That is, it displays a fractional period with respect to 
the unit, which is the Heisenberg time. This period has the form T = 1/c, with c an 
integer whose specific value is linked with the value of the rank of the interaction k of 
the residual interaction. 

These fractional periodic revivals are observed when some specific matrix elements 
of the fc-body interaction matrix v r>s dominate over the rest. In fact, these matrix ele- 
ments correspond to eliminating precisely c particles from one single-particle state, and 
creating them in the other one. By definition of the £-body interaction, c is bounded as 
1 < c < k. Therefore, for k = 2 we may move at most two particles from one single- 
particle state to the other one. When the associated matrix elements somehow dominate 
the perturbation, fractional periodic revivals are observed with period 1/2 in units of 
tH (see Fig. 3). Likewise, for k = 3 there are perturbing terms that involve moving up 
to one, two or three particles. The relative dominance of these terms is related to the 
observation of periodic revivals with periods 1, 1/2 or 1/3 in Heisenberg-time units, 
respectively (see Fig. 4). We note that these results do not depend on the time-reversal 
invariance of the ensemble, given by /3 . 

The results discussed above are illustrated in Figs. 3 and 4. In Fig. 3 we present 
the fidelity decay for short times (in units of tn) of two realizations of the ensemble 
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FIGURE 3. Fidelity decay for small times of two specific realizations of the ensemble forfc = 2 (« = 128 
and X = 10~ 6 ). During the freeze of fidelity, periodic revivals are displayed, with period 1 (dotted curve), 
and of period 1/2 (solid curve), in units of the Heisenberg time. The inset presents the same results for 
longer times. 

for k = 2 in = 128 and X = 1(T 6 ). The dotted curve displays a case comparable 
with the ensemble-average fidelity, where the terms moving one particle dominate the 
perturbation; the solid curve illustrates a case where the oscillations of period 1/2 

a offdl£l2 

appear. In the latter case, the dominating two-body interaction H k corresponds to 

the term (d\) 2 (d2) 2 + (d^) 2 (di) 2 . Figure 4 illustrates the k = 3 case. As discussed above, 
in the case of three-body interactions it is possible to observe oscillations of period 1 
(dashed line), 1/2 (dotted line) or 1/3 (solid line), depending on the specific £-body 
matrix elements that dominate the interaction. 

4. CONCLUSIONS 

In this paper we have presented numerical results on the fidelity decay of the two- 
level bosonic fc-body embedded ensembles of random matrices, for k = 2 and k = 
3, considering the ensemble-averaged fidelity and comparing it with some individual 
realizations of the ensemble. The reference Hamiltonian is diagonal in the occupation 
number basis, where the residual interaction is traceless. We have obtained the relevant 
scaling laws, and observed the fidelity freeze, which displays oscillations of period 1 in 
units of the Heisenberg time. 

We have also presented preliminary numerical results on the existence of fractional 
periodic revivals (in units of the Heisenberg time) during the freeze of fidelity, for some 
specific members of this ensemble. Expressed in units of the Heisenberg time, the period 
of these revivals is T = 1 /c, where c is an integer number related to the specific k- 
body interactions which dominate the perturbation. The occurrence of these fractional 
periodic revivals is somewhat rare with respect to ensemble average; we conjecture 




FIGURE 4. Same as Fig. 3 for k = 3 (n = 128 and A = 1(T 6 ), displaying cases where the periodicity of 
the revivals is 1 (dashed curve), 1 /2 (dotted curve) or 1 /3 (solid curve), in units of the Heisenberg time. 
The inset shows the same results for longer times. 



that this is related to the relative number of terms in the perturbation which move c 
particles from one level to the other. Our results indicate that the appearance of the 
fractional periodic revivals is independent of the time-reversal invariance. These results 
may be interesting for the understanding — and even the measurement — of three-body 
interactions in two-component Bose-Einstein condensates [15]. 
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